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Abstract 

Near-surface soil moisture is a critical component of land surface energy and water balance 
studies encompassing a wide range of disciplines. However, the processes of infiltration, runoff, and 
evapotranspiration in the vadose zone of the soil are not easy to quantify or predict because of the 
difficulty in accurately representing soil texture and hydraulic properties in land surface models. 
This study approaches the problem of parameterizing soils from a unique perspective based on 
components originally developed for operational estimation of soil moisture for mobility 
assessments. Estimates of near-surface soil moisture derived from passive (L-band) microwave 
remote sensing were acquired on six dates during the Monsoon ’90 experiment in southeastern 
Arizona, and used to calibrate hydraulic properties in an offline land surface model and infer 
information on the soil conditions of the region. Specifically, a robust parameter estimation tool 
(PEST) was used to calibrate the Noah land surface model and run at very high spatial resolution 
across the Walnut Gulch Experimental Watershed. Errors in simulated versus observed soil 
moisture were minimized by adjusting the soil texture, which in turn controls the hydraulic 
properties through the use of pedotransfer functions. By estimating a continuous range of widely 
applicable soil properties such as sand, silt, and clay percentages rather than applying rigid soil 
texture classes, lookup tables, or large parameter sets as in previous studies, the physical accuracy 
and consistency of the resulting soils could then be assessed. 

In addition, the sensitivity of this calibration method to the number and timing of microwave 
retrievals is determined in relation to the temporal patterns in precipitation and soil drying. The 
resultant soil properties were applied to an extended time period demonstrating the improvement in 
simulated soil moisture over that using default or county-level soil parameters. The methodology is 
also applied to an independent case at Walnut Gulch using a new soil moisture product from active 
(C-band) radar imagery with much lower spatial and temporal resolution. Overall, results 
demonstrate the potential to gain physically meaningful soils information using simple parameter 
estimation with few but appropriately timed remote sensing retrievals. 
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Abstract 


Near-surface soil moisture is a critical component of land surface energy and water 
balance studies encompassing a wide range of disciplines. However, the processes of 
infiltration, runoff, and evapotranspiration in the vadose zone of the soil are not easy to 
quantify or predict because of the difficulty in accurately representing soil texture and 
hydraulic properties in land surface models. This study approaches the problem of 
parameterizing soils from a unique perspective based on components originally developed 
for operational estimation of soil moisture for mobility assessments. Estimates of near- 
surface soil moisture derived from passive (L-band) microwave remote sensing were 
acquired on six dates during the Monsoon ’90 experiment in southeastern Arizona, and used 
to calibrate hydraulic properties in an offline land surface model and infer information on the 
soil conditions of the region. Specifically, a robust parameter estimation tool (PEST) was 
used to calibrate the Noah land surface model and run at very high spatial resolution across 
the Walnut Gulch Experimental Watershed. Errors in simulated versus observed soil 
moisture were minimized by adjusting the soil texture, which in turn controls the hydraulic 
properties through the use of pedotransfer functions. By estimating a continuous range of 
widely applicable soil properties such as sand, silt, and clay percentages rather than applying 
rigid soil texture classes, lookup tables, or large parameter sets as in previous studies, the 
physical accuracy and consistency of the resulting soils could then be assessed. 

In addition, the sensitivity of this calibration method to the number and timing of 
microwave retrievals is determined in relation to the temporal patterns in precipitation and 
soil drying. The resultant soil properties were applied to an extended time period 
demonstrating the improvement in simulated soil moisture over that using default or county- 



level sell parameters. The methodology is also applied to an independent ease at Walnut 
Guleh using a new soil moisture product from active (C-hand) radar imagery with much 
lower spatial and temporal resolution. Overall, results demonstrate the potential to gain 
physically meaningful soils information using simple parameter estimation with few but 
appropriately timed remote sensing retrievals. 
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1. Introduction 


Soil moisture remains an essential yet elusive component of Earth system seienee research 
across a wide range of scales and applications. In addition to impacting agriculture, water 
resource management, and extreme events such as flooding and drought, the day-to-day 
variability in soil moisture on field to global scales is an important quantity for atmospheric 
modeling and prediction. In fact, the accuracies of climate, mesoseale, boundary layer, land 
surface, and hydrologic models are ultimately dependent on proper treatment and simulation of 
the state and transfer of water and heat at the land surface (Koster et al. 2004; Findell and Eltahir 
2003; Berbery et al. 2003; Betts et al. 2003; Betts 2000). 

Unfortunately, soil moisture is not as easily measured or observed as atmospheric properties 
such as temperature, humidity, and wind speed. For example, in situ or remotely-sensed 
observations of soil moisture for initialization, update, and validation purposes are not yet 
available on the scales of most models. Observations are generally confined to short-term field 
experiments, many of which have highlighted the heterogeneous nature of soils in terms of water 
content and texture (Mohanty et al. 2002). Indirect estimates of soil moisture can be obtained 
using thermal infrared measurements (Carlson et al. 1995), but require a priori information on 
the surface characteristics. As an alternative, passive and active microwave remote sensing 
methods have had the greatest success in estimating soil moisture in a temporally and spatially 
consistent manner (Thoma et al. 2006; Moran et al. 2004; Hollenbeck et al. 1996). 

Recent studies have noted that the most successful and promising approach to estimating 
soil moisture continuously over time and space must include a combination of remote sensing 
and modeling (Entekhabi et al. 1999; Houser et al. 1998). The majority of land surface models 
(LSMs) require soil hydraulic parameters to solve for the transport of moisture within the soil 
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using Richards’ (1931) formulations. These parameters are often derived from soil texture 
information, but due to the heterogeneous nature of soils and lack of detailed maps of soil 
properties, soil parameterization schemes are often crude, inflexible, or inappropriate. Further, 
LSMs have been shown to potentially be more sensitive to the choice of soil hydraulic properties 
or soil texture data than to atmospheric forcing or surface characteristics (Gutmann and Small 
2005; Pitman 2003). 

Because of these difficulties, numerous attempts have been made to optimize 
LSM parameters using observations of state variables such as soil moisture and surface 
temperature as constraints (Hogue et al. 2005; Liu et al. 2004; Hess 2001; Gupta et al. 1999 ). 
While these studies highlight the potential for parameter estimation techniques to derive large 
sets of ‘effective’ parameters and diagnose specific model weaknesses, little has been gained in 
terms of acquiring physieally-meaningful or hydraulically consistent estimates of individual 
parameters. Because of the complexity and number of estimation techniques and parameter sets 
employed in these studies, it remains difficult to infer or derive any parameter information that 
could be applied to other independent studies or models. 

With these issues in mind, this paper examines the potential use of passive and active 
microwave retrievals of near-surface soil moisture to calibrate a LSM and infer a physically- 
meaningful and consistent set of soil hydraulic parameters, using a combination of high- 
resolution land surface modeling and parameter estimation. The experimental design of this 
work was originally developed for the purpose of estimating troop and vehicle mobility for the 
United States Army based on operational soil moisture prediction from a very limited set of input 
data (Army Remote Moisture System; ARMS; Tisehler et al. 2006). Here, we have tested and 
extended ARMS to assess the ability of parameter estimation techniques to minimize inherent 
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model error, yet still provide information on difficult to obtain soil properties over tbe Walnut 
Guleh Environmental Watershed (WGEW) in Arizona. 

Accordingly, Section 2 summarizes tbe current state of knowledge of the many components 
of the ARMS project including soil parameterizations in LSMs, microwave remote sensing of 
soil moisture, and parameter estimation. In Section 3, the models, site, and remote sensing data 
employed in this study are described. Results of the calibration experiments are presented in 
Section 4, including an evaluation of the optimized parameters and sensitivity to temporal 
sampling of remote sensing. Finally, Section 5 discusses the limitations and applicability of the 
results, including suggestions for the future utility of physically meaningful parameters in LSMs. 

2. Background 

a. Soil Parameterizations in LSMs 

The in f lu ence of near-surface soil moisture on the partitioning of surface turbulent fluxes 
from offline LSMs to folly coupled global climate models has been well-documented (e.g., 
Braun and Sehadler 2005; Ek and Holtslag 2003; Santanello and Carlson 2001; Cuenca et al. 
199b; Sun and Bosilovieb 199b; Ek and Cuenca 1994; Jacobs and DeBruin 1992). In order to 
simulate the evolution of moisture in the soil, a set of soil hydraulic parameters are combined 
with expressions (known as soil moisture characteristic curves) relating soil moisture (0 with 
matric potential ( i//), and soil moisture with hydraulic conductivity (K). The expressions derived 
by Brooks and Corey (1964) and Campbell (1974) are most commonly used in meteorological 
coupled models, while the van Genuehten (1980) functions based on a different set of soil 
measurements are used for more detailed soil and hydrological models. A full description and 
evaluation of these functions can be found in Braun and Sehadler (2005). 
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The three forms of the characteristic corves above depend on a set of 4 (Campbell, 1974} or 
5 (Brooks and Corey, 1964; van Genuehten, 198Q) hydraulic parameters, which are a function of 
the soil composition and structure. These parameters include the saturated matric potential ( i// s ; 
aka “bubbling” or “air entry”), the saturated hydraulic conductivity (K s ), the saturated soil 
moisture content (porosity; 60, the residual soil moisture content (60, and the pore size 
distribution index (b). Unfortunately, estimating these parameters consistently and accurately 
has proven difficult even for identical soils measured under controlled laboratory conditions. 
Further, studies have shown that LSM simulation of soil moisture can be more dependent upon 
the specification of hydraulic parameters than atmospheric forcing or surface conditions 
(Gutmann and Small 2005; Santanelloand Carlson 2001). 

To mitigate these differences and acquire a somewhat standard set of parameters for LSM 
applications, ‘hulk’ parameters have been derived that are based on soil type. The results of 
Clapp and Homberger (CH; 1978), Rawls et al. (1982), and Coshy (1984) and provide the most 
extensive and commonly employed lookup tables of hydraulic parameters for LSMs, with 
atmospheric-based applications favoring CH and Cosby and soil hydrology models employing 
the Rawls parameters. Unfortunately, parameter lookup tables are only as accurate as the 
available soil texture type information and provide an “average” value of each parameter for 
each soil type. High-resolution soil texture maps remain difficult to obtain, particularly on 
regional and global scales, and there is little flexibility between soil types or for mixed soils 
despite that larger differences in soil properties have been observed within a certain soil type 
than between types (Gutmann and Small 2005; Soet and Strieker 2005; Feddes et al. 1995). 

To bridge the gap between rigid soil textural classes and the heterogeneous nature of soils, 
numerous pedotransfer functions (PTFs) have been developed (Sobieraj et al. 2001). The most 
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commonly used ‘class’ PTFs relate discrete soil types to hydraulic parameters and are the basis 
upon which lookup tables are used in L&M and meteorological modeling applications. 
‘Continuous’ PTFs are more detailed and relate measurable soil properties such as percent of 
sand and clay, porosity, and bulk density to hydraulic properties using regression equations 
derived from soil samples. These functions are continuous without bounds, and therefore allow 
more flexibility and independence in parameter values than those from lookup tables. More 
importantly, continuous PTFs that are able to reproduce areal averaged conditions in LSMs have 
been shown to scale linearly in space and therefore could be used to infer spatially-aggregated 
hydraulic parameters. Although the advantages of continuous over class PTFs has been 
demonstrated for hydrologic models (Soet and Strieker 2003), continuous PTFs are not routinely 
employed in LSMs (except for CLM, give reference) or atmospheric models where the broad 
definition and application of soil types still dominate the simulation of soil moisture. 
b. Parameter Estimation 

An alternative to specifying highly-uneertain soil hydraulic parameters in LSMs is to use 
parameter estimation and model calibration techniques. For example, a relatively simple and 
well-established parameter estimation model (PEST; Doherty 2004) has been used by a number 
of scientific disciplines to optimize model parameters given limited observations of fundamental 
output variables. For example, by adjusting soil porosity in a LSM until the difference in 
simulated versus observed soil moisture is minimized (through a specified objective function), an 
LSM can be calibrated using PEST. 

In recent years, more sophisticated estimation techniques have been developed to estimate 
large and diverse sets of parameters. Liu et al. (2003) used a multi-objective technique for 
offline and partially-coupled LSMs to examine the pathways by which a deficiency in the model 
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physics impacts coupled and uncoupled simulations. Following this work, Liu et al. (2004) 
performed controlled parameter estimation studies of offline and partially-coupled models and 
examined the effects of including atmospheric (in addition to soil and vegetation) parameters in 
the optimization. Hogue et al. (2005) investigated the transferability of large optimized 
parameter sets in an offline LSM across varying surface conditions and time periods, and 
concluded that parameter optimization needs to be site-specific for best results, and should be 
recalibrated for changes in seasons or over longer time intervals. 

Scott et al. (2000) performed soil hydraulic parameter estimation using the Hydrus soil 
moisture model at two sites in die Walnut Gulch Environmental Watershed (WGEW) in Arizona. 
While the focus was on the vertical distribution of soil moisture and recharge at these points 
alone, their results show that the model was least sensitive to K sat and most sensitive to porosity 
and b, which is consistent wife other studies. Scott et al. (2000) also stress feat fee derived 
parameters are ‘effective’ in nature, compensating for errors in the soil physics of the model, and 
feat further research is needed to assess the limitations of parameter estimation across spatially 
heterogeneous and distributed watersheds. 

Overall, parameter estimation studies have focused on large sets of parameters and complex 
algorithms that require a great deal of computational time. From these studies, it could also be 
argued that the bulk of the work done to this point has been focused on ‘model calibration’ rather 
than parameter estimation, particularly when there is significant model error accounted for in the 
optimized parameters. It is important to note feat the research presented here differs from such 
multi-objectives techniques, and is focused solely on calibrating a physically meaningful set of 
soil hydraulic properties that improve soil moisture simulated by a LSM. 

c. Remote Sensing of Soil Moisture 
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Due to the limited nature of available soil instrumentation and measurement techniques 
(e.g., theta probe, TDR, Vitel probe, gravimetric), a spatially continuous and reliable network of 
soil moisture measurements that could be used to initialize and evaluate LSMs does not exist. 
As a result, passive microwave (L-band; 1.4 GHz) estimation of soil moisture has been has been 
explored a great deal using instruments such as NASA’s push broom microwave radiometer 
(PBMR; Sehmugge et al. 1988). Due to the high spatial resolution required at this wavelength 
(21 cm), passive microwave radiometers are typically flown on aircraft where they have shown a 
great deal of promise in estimating soil moisture across varying surface conditions (Mattikalli et 
al. 1998; Burke et al. 1997; Hollenbeck et al. 1996). Changes in the dielectric constant the top 5 
cm of soil are due to changes in the relative water content, and are evident in the brightness 
temperature measured by the sensor. 

More recently, techniques have been developed to estimate soil moisture using active 
microwave remote sensing (C-band; 5.3 GHz). Because of the shorter wavelength (5.6 cm), 
active sensors can placed aboard satellite platforms and potentially acquire high resolution 
estimates of soil moisture when combined with empirical and physical models (Thoma et al. 
2006). To date, there have been mixed results using radar remote sensing to estimate soil 
moisture due to the sensitivity of low frequency baekseatter to the nature and degree of surface 
interactions and, consequently, the degree of signal correction required (see also review by 
Moran et al. 2004). 

Recently, Thoma et al. (2006) have developed an image differencing technique for active 
remote sensing that shows promise in eliminating much of the noise in C-band radar data. This 
‘delta index’ method requires a single reference (dry) image to compare with a separate (wet) 
image over the same domain (assuming no other changes in surface characteristics), thereby 
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isolating the change in baekseatter due to soil moisture variations. This method acts to minimize 
errors due to surface roughness effects using filtering techniques to reduce the amount of speckle 
that is common in radar imagery (particularly in regions of high rock fragment). 

The delta index is defined as follows. 


delta index = 



( 1 ) 


where o&y is the baekseatter (db) from a dry radar image, and o wet is the radar baekseatter (db) 
from the identical pixel location in a wet image. The delta index has been shown to have a near 
linear (1:1) relationship with volumetric soil moisture, and is particularly applicable to semi-arid 
regions where a spatially-uniform dry reference image can be acquired (Thoma et al. 2006). 

d. Estimation of Soil Hydraulic Properties 

Since the development of L-band passive microwave soil moisture retrievals, numerous 
studies have attempted to use a combination of remote sensing imagery, LSMs, radiative transfer 
(emission) models, and observations to infer soil hydraulic properties. For example, van de 
Griend and O’Neill (1986) demonstrated that independent measurements of soil moisture from 
microwave remote sensing and the thermal inertia of the soil can be related to hydrologic 
properties of loamy sand soils during an 11 -day dry down period. This work was extended by 
Canrillo et al. (1986) using a combination of models and measurements for three distinct soil 
types under highly controlled plot-scale conditions. They calibrated a soil model (hydraulic 
properties) until a coupled microwave emission model best matched the observations of L-band 
microwave brightness temperature over a 3-day drydown. Overall, Camillo et al. (1986) suggest 
that a wider range of soil moisture conditions than those observed here may improve results by 
better capturing the functional drying curves represented by the soil model parameterizations. 
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Following the work of Camilla et al. (1986), Burke et al. (1997) and (1998) used a coupled 
land surface-microwave emission model in conjunction with radiometer measurements to infer 
soil properties for hare and vegetated soil plots. Performed over a IQ-day period with a primarily 
sandy soil and hare soil, com, and soybean canopies, hydraulic parameters were adjusted to 
match the emission model output with L-band radiometer measurements. In agreement with 
other studies, the model was found to be least sensitive to K sat and most sensitive to b, and for the 
com and soybean plots vegetation parameters such as leaf area index and root density were 
significant. Overall, these laboratory and point-scale studies point towards the future use of 
PTFs rather than a one-at-a-time parameter estimation approach to acquire spatially-distributed 
soil properties over watersheds, and suggest that an intensive period of microwave images be 
acquired to capture significant soil drydowns. 

Feddes et al. (1993) examined the use of microwave measurements of soil moisture, 
temperature, and albedo to calibrate and infer soil hydraulic properties. They found that the 
‘effective’ soil parameters for the LSM could be derived using this approach. However, their 
method also required a great deal of measurements and parameters, such as evaporation and in- 
situ soil moisture at multiple depths, thereby limiting its application to highly controlled and 
plot-scale experiments. 

Hollenbeck et al. (1996) used PBMR estimates of soil moisture to infer soil conditions 
during large-scale field experiment at HAPEX-Sahel. Two PBMR brightness temperature 
images, two days apart, were used to calculate the relative change in soil moisture following a 
precipitation event to infer soil hydraulic properties over the watershed. Though their results 
were entirely based on qualitative image-differencing, they were able to isolate the impact of soil 
properties on the image differences versus that of initial soil moisture, land cover, and rainfall 
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distribution. The impact of antecedent precipitation is suggested for further study, as it greatly 
impacts the stage of soil drying being monitored by the PBMR. 

Finally, Mattikalli et al. (1998) tested the laboratory results of Ahuja (1993), who 
demonstrated that K sat could be derived using remotely sensed estimates of 2-day changes in soil 
moisture. They concentrated on calibrations of a hydrologic model for 3 layers of soil moisture 
and parameters across 13 sites in the Little Washita, OK watershed. A significant qualitative 
correlation between spatial maps of brightness temperature, soil moisture, and soil texture give 
validity to the strong relationships between microwave measurements and soil type and 
properties for this region. Although treatment of the remaining hydraulic parameters, spatial 
distribution of K sah or a detailed evaluation against typical K sa , values was not presented, this 
study confirms the theoretical framework by which a more comprehensive approach to 
estimating these parameters can be based. 

e. Summary 

These studies have demonstrated that the strong link between microwave remote sensing 
and soil moisture (that is ultimately controlled by hydraulic parameters) can provide a pathway 
to improve LSM soil physics and parameterizations. While these works have provided a strong 
physical and methodological foundation by which to address these issues, each has limitations in 
te rms of scope and applicability that can now be improved upon by taking the suggested next 
steps and utilizing new approaches and data. Specifically, this paper will bridge the gaps 
between and extend previous studies by: 

1} Determining the ability of parameter estimation to calibrate a LSM and to infer physically 
meaningful estimates of soil hydraulic properties using pedotransfer functions and 
microwave remote sensing of soil moisture at high spatial and temporal resolution; 
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2) Testing the sensitivity of the calibration process and retrieved properties on precipitation 
and soil drydown patterns using temporal sampling of remote sensing imagery; and 

3) Applying the retrieved soil parameters to an independent dataset, and assess the ability of 
a new image differencing technique of estimating soil moisture from active microwave 
remote sensing to he used in the calibration process. 

3. Methodology and Data 

a. ARMS Background 

The Army Remote Moisture System (ARMS; Tisehler et al. 2006) project is an ongoing 
collaboration between the U. S. Army Corps of Engineers, U. S. Department of Agriculture, 
NASA-GSFC, and the University of Wyoming. The goal of this work is to provide improved 
operational estimates of soil moisture and hydraulic properties as inputs to decision-making 
models based on factors such as troop and vehicle mobility and landing strip suitability. The 
three main components of ARMS are 1) high-resolution microwave remote sensing of soil 
moisture, used to calibrate a 2) land surface model by optimizing hydraulic properties through 3} 
parameter estimation. The ultimate goal of ARMS is to be able to use limited site information 
and radar-based soil moisture retrievals to calibrate an LSM for any location in the world and 
enable soil moisture and properties to be more accurately simulated going forward. While this 
study is focused on a semiarid testbed in Arizona, ARMS is also being tested at other diverse 
locations across the U. S. (OK, GA, and CO). 

b. Site Information 

The Walnut Gulch Experimental Watershed (WGEW) is a located in southeastern Arizona, 
covering 148 km 2 of semi-arid grassland and shrub covered rangeland. The detailed 
instrumentation and long record length of the datasets available in this region have made the 
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WGEW the focus of many hydrological, meteorological, and remote sensing studies. Most 
notably, the Monsoon ’90 field experiment (M9Q; Kustas et al. 1991} was conducted in this 
region in July and August of 1990, and included the deployment of eight Metflux sites across the 
watershed that measured standard meteorological data as well as land cover, soil moisture, and 
soil property information. (Fig. 1} 

Overall, the conditions throughout the WGEW are dominated by the summer monsoon 
season of July and August, when the bulk of the annual 250-500 mm rainfall (mainly convective) 
occurs. Rainfall events during the monsoon period are typically <10 mm and only influence the 
top 10 cm of soil before being quickly returned to the atmosphere through ET within 3 days 
(Kure and Small 2004). This means that the near-surface soil moisture is the only variable 
reservoir of moisture in this region. During the period from April-July, the soils often reach a 
desiccated state before the onset of the monsoonal precipitation. Land cover consists mainly of 
open shrubland (< 30% cover) in the western half of the WGEW, and grass cover (< 50% cover) 
in the east. 

At each Metflux site, standard meteorological variables were measured at 20-minute 
intervals. Precipitation measurements were derived from a dense 98-gauge network covering the 
entire watershed, from which spatially interpolated rainfall estimates (useful for modeling 
applications) have been generated using a variety of techniques (Garcia et al. 2006; Houser et al. 
1998). Two supersites were furnished with additional instrumentation: Lucky Hills (LH) 
located in the shrub dominated north-central part of the domain, and Kendall located in the 
grasslands of the east. Soil moisture in the upper 5 cm layer was estimated at each site using 
multiple gravimetric measurements, and vertical profiles of soil moisture were estimated at 
Kendall and LH using time domain refleetrometry, giving estimates from 5-50 cm in depth. 
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In addition to the M90 period, extensive meteorological and flux measurements at the LH 
and Kendall sites and eo-located precipitation and soil moisture measurements at the rain gauge 
stations have been acquired for the 2002-2004 period. This dataset is temporally more extensive 
than that of M90 and provides an independent dataset to evaluate LSM simulations as well as 
new remote sensing data and parameter estimation techniques. 
c. Remote Sensing of Soil Moisture 

Passive microwave remote sensing measurements (L band; 21 -cm) of brightness 
temperature were made over a significant portion of the WGEW during the M90 experiment 
using NASA’s push broom microwave radiometer (PBMR). As mentioned, the PBMR has been 
used extensively for measuring soil moisture across arid and semi-arid regions. Flights taken 
during the Monsoon ’90 experiment are described in detail in Sehmugge et al. (M9Q; 1994), and 
during the Hydrologic Atmospheric Pilot Experiment in 1992 (HAPEX-Sahel; Hollenbeck et al. 
1996). From this dataset, six daily estimates of near-surface soil moisture are available both 
before (D0Y 212) and after (DGY 214, 216, 217, 220, and 221) the onset of precipitation. The 
PBMR data was resampled to 40 m resolution and mapped to a UTM grid that covers a subset of 
the WGEW that includes all 8 Metflux sites. Sehmugge et al. (1994) showed that brightness 
temperature measurements correlated well with both rainfall and 0-5 cm soil moisture measured 
at the sites. Figure 2 shows the PBMR and gravimetric estimates of soil moisture at the Kendall 
and LH sites along with gauge-interpolated precipitation during the M90 period. These plots 
highlight the desiccated soil conditions before the first and most intense rainfall event on DOY 
2 IS, and more importantly how the PBMR images captures the period of rapid soil drying 
thereafter. 
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For the 2002-4 period active, or radar, microwave remote sensing (C band; 5.6cm) 
measurements were acquired from RADARSAT-1 imagery. A reference (dry) image was taken 
on 19 January 2003 and combined with images during the monsoon period (29 July, 22 August, 
and 15 September 2003) to derive soil moisture using the delta index approach (Eqn. 1). The 
three resulting 0-5 cm soil moisture estimates cover a 6 week period spanning an extended 
(seasonal) drydown period immediately following rainfall. The nominal resolution of 
RADARSAT-1 is 7 meters and covers the entire WGEW domain, but the raw baekseatter data 
was further processed and filtered as discussed by Thoma et al. (2006) to 210 and 280 meter 
resolution to reduce the effects of speckle. 
d.LSM 

The Noah land surface model (Chen et al. 1996) was originally developed from the land 
component of the Oregon State University 1-D planetary boundary layer model (OSU; Mahrt 
and Ek 1984), and is currently employed as the land surface scheme in NCEP’s operational 
version of the Weather Research and Forecasting nonhydrostatic Mesoseale Model (WRF- 
NMM). Extensive evaluations and discussion of the Noah physics and comparisons to other 
LSMs has been performed by Roboek et al. (2003) among others. The hydrology within the 
Noah model is handled by a Richard’s equation formulation governed by the Campbell functions 
(1974). Hydraulic parameters are typically estimated from prescribed soil types based on the 
lookup tables and regression results of Cosby et al. (1984). Noah is typically used across a wide 
range of scales offline and for meso- to global-scale meteorological applications when coupled 
with atmospheric models. 

The offline version (2.6) of Noah was configured to run at 40 m resolution over the WGEW. 
Forcing data was acquired from the LH site and applied uniformly to each pixel within the 


15 



domain including downward shortwave and longwave radiation, air temperature, specific 
humidity, and wind speed at a specified reference height (6 m), and surface pressure. As- 
suggested by Houser et al. (1998), the impacts of using of a single forcing dataset applied to the 
entire watershed are minimal so long as spatially distributed precipitation is accounted for, which 
is the ease here. Indeed, si m ulations were performed using the Kendall forcing data in place of 
LH, and resulted in changes in surfaces fluxes that were less than the instrument error. To 
ensure that the most detailed precipitation data was included, rainfall from 84 of the 98-gauge 
network was broken down and interpolated in time and space at 20-minute and 40-meter 
resolution. Garcia et al. (2006) provide a detailed description and evaluation of two methods of 
interpolating raingauge data over the WGEW, and for the large number of gauges available here 
both the inverse distance weighting (IDW) an d multi-quadric biharmonic (MQ-B) schemes work 
equally well for precipitation interpolation. 

In Noah, vegetation parameters for each grid cell are derived from land cover maps (using 
lookup tables). In a similar manner, ten soil parameters are typically derived at each grid cell 
from soil texture maps from the Food and Agricultural Organization of the United Nations 
(FAO; FAO-UNESCO 1984), State Soil Geographic Database (STATSGO; USDA 1994), or 
Soil Survey Geographic Database (SSURGO; USDA 2002) datasets. Moeko et al. (2006) 
performed a thorough evaluation of the impacts of varying inputs of land cover, soil type, and 
precipitation on soil moisture simulations. In this study, we use the best available land cover 
data (UMD; Hansen et al. 2000), elimatok>gieally-derived values of albedo and vegetation 
fraction, and MQ-B precipitation forcing as input to the Noah model. 

A standard 4-layer soil profile was used in Noah with a top layer of 5 cm that matched the 
representative depth of in-situ and remote sensing soil moisture measurements. The sensitivity 
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of Noah to varying numbers (up to 20) and depths of sod layers was examined in detail, and 
results showed that adding additional layers and adjustment of vegetation rooting depth did not 
significantly alter or improve results for this semi-arid region. Also, after careful calibration 
with observed soil moisture values in the WGEW and the results of previous studies (Scott et al. 
2000), the prescribed minimum value of soil moisture in Noah (i.e. wilting point) was lowered 
from 0.05 to 0.02 m 3 m' 3 . 

In an effort to ensure consistency and add flexibility within soil types and hydraulic 
parameters, PTFs were incorporated into Noah for this study. Specifically, the PTFs derived by 
Cosby et al. (1984) require only percentages of sand and clay to derive the hydraulic parameters, 
which are then independent of soil texture classes or averaging. As discussed, previous studies 
have pointed to the potential advantages of PTF approaches over discrete soil types and lookup 
tables in LSMs. Though based on the identical soil samples and data of the default Noah lookup 
tables, these PTFs ensure that a fell range of soil parameter values based on soil composition are 
derived in a realistic and consistent maimer. 
e. Parameter Estimation 

The Parameter Estimation model (PEST; Doherty 2004) is a widely-used tool for examining 
sensitivities and estimating parameters in models spanning a wide range of applications. In 
particular, the ability of PEST as a model-independent estimation technique to link with any type 
of LSM using flexible parameter, observation, and convergence criteria make it optimal for use 
in this study. Here, PEST was configured to run as a parent model to Noah (PEST-Noah), where 
it evaluates and minimizes an objective function based on the differences between simulated and 
observed soil moisture as follows 

RMSE = 
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where N 0 bs is the number of PBMR observations used in the calibration and Qwoak and &x,pbmr are 
Noah simulated and PBMR observed 0-5 cm soil moisture at each observation time. Until the 
convergence criteria are met, PEST iterates and adjusts the Noah parameters (sand, silt, and clay 
percentages), evaluates if this decreases the model error, and adjusts the parameters accordingly. 
Extensive testing of PEST-Noah has shown that there can be on the order of 2-20 optimizations 
requiring up to 200 model runs before PEST converges in some eases, depending on how far 
apart the initial parameters are from their optimal values. 

To ensure the accuracy and repeatability of PEST-Noah simulations, synthetic twin 
experiments were conducted. A single (control ease) Noah simulation of the M90 period with 
prescribed sand, silt, and clay values was performed, from which the 0-5 cm soil moisture output 
was extracted on six dates corresponding to the PBMR observations. PEST-Noah was then run 
using Noah soil moisture output (in place of PBMR) as observations until the optimized sand, 
silt, and clay percentages were found. This process was repeated at the Kendall and LH sites 
using vastly different soil textures and moisture contents, and in each case PEST-Noah returned 
the precise sand, silt, and clay values prescribed in the control case. This gives confidence that 
PEST-Noah can be run for a variety of conditions at WGEW and that the results are unique. 

4. Results 

a. M90 Calibration Experiments 

Simulations were performed during the M90 period from 23 July - 9 August 1990 that 
encompasses the 6 PBMR overpasses. This period allows ample time for the model to 
equilibrate to the very dry initial conditions leading up to the first PBMR image (31 July), and 
before the onset of the monsoon and the first significant precipitation event of the season (2 
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August}. The model was run using a 20-minute timestep, and output was generated every six 
hours. 

i. Metflux Sites 

PEST-Noah was run at each of the 8 Metflux sites using the closest 40m PBMR pixel to 
each site on the 6 observation dates. Figure 3 shows the simulated soil moisture at the Kendall 
and LH sites before (FAO soils) and after (PEST) calibration of sand, silt, and clay along with 
corresponding PBMR and in-situ gravimetric measurements. Despite the differences in 
magnitude and drydown patterns exhibited between the sites, PEST is able to fit the simulated 
soil moisture to the observations. Also evident is the significant improvement in simulations 
using calibrated soil properties compared with those from default lookup table approach. 

The RMSE and bias in simulated versus observed (PBMR) soil moisture for all 8 Metflux 
sites are shown in Figure 4. The ARMS requirement of 5 percent (volumetric) accuracy in soil 
moisture prediction is easily satisfied at all the Metflux sites when using the PEST-Noah 
calibration, with overall RMSE and bias values less than 2 percent. In particular, the bias in the 
default Noah simulation using FAO soils has been greatly reduced using PEST to near zero for 
most locations. Examination of each individual site’s improvement in simulated soil moisture 
(similar to that shown in Figure 3) makes it clear that PEST-Noah primarily acts to reduce the 
bias by adjusting the overall magnitude and dynamic range (using soil texture) to match 
observations. 

Given the accuracy of the calibrated soil moisture at each site, it is useful to assess the 
potential utility and accuracy of optimized soil textures as well. Figure 5 shows a comparison of 
the optimized sand, silt, and clay percentages at each site versus those measured during the M90 
experiment by Schmugge et al. (1994). The optimized soil textures suggest a primarily sandy 
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soil, and are similar to the observed soils with greater than 65 percent sand and less than 10 
percent clay content. Optimized values of silt content are lower than those observed, but this is 
likely due to the setup of PEST-Noah where silt is actually a dummy variable and, more 
importantly, that the PTFs are only a function of sand and clay content. Therefore, adjusting silt 
in the simulations had zero impact on the soil moisture and value of the objective function 
evaluated by PEST. 

It is also important to consider the range and magnitude of hydraulic properties resulting 
from the different soil textures. Table la presents the optimized values of sand, silt, and clay for 
the Metflux sites and corresponding hydraulic properties derived from the PTFs in Noah. For 
comparison, Table lb lists observed soil textures from Sehmugge et al. (1994) and the North 
American Monsoon Experiment (NAME; Higgins et al. 2006) (sites 1 and 5 only), and hydraulic 
properties estimated using the Noah PTFs. Overall, there is relatively little variation in 
properties estimated from PEST-Noah across the sites despite variation in sand content (73-100 
percent). Similarly, hydraulic properties derived from observed textures exhibit a small range, 
although their magnitude differs slightly from the PEST-Noah values due to the lower sand 
percentage (66-80 percent). 

PEST-Noah suggests a slightly more sandy soil over the region than observed, but the 
differences in the parameters that control the soil moisture dynamics are not as significant as the 
soil textures might indicate. Plausible values are estimated by PEST-Noah for each parameter 
when compared with observations (even at site 7 where an estimate of 100 percent sand may not 
appear physically realistic, but the resulting properties are within measured values). This is due 
to the slowly varying relationships between soil texture and hydraulic properties for high sand 
contents governed by the PTFs. 
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Evidence that the spatial differences in soil texture from PEST-Noah also have physical 
meaning is present in Fig. 5 and Table 1. For example, sites 2 and 7 are the extremes in 
optimized soil textures from lowest (72 percent) to highest (100 percent) sand percentage. 
Measurements at these sites also support that site 7 is sandiest, and site 2 has the highest mixture 
of silt and clay. PEST-Noah follows the same trend at site 2, estimating a significant silt 
percentage which, in effect, allows the sand percentage to be lower. Inspection of the soil drying 
patterns observed at these sites during M90 confirms that at site 2 there is higher overall soil 
moisture and a more consistent but slow drydown compared with site 7. These characteristics 
and subtle differences in the drying curves are responsible for PEST-Noah estimating a less 
sandy soil at site 2, which is representative of the observed local soils. 

To get a better feel for the physical applicability of the parameters themselves. Table 2 
presents hydraulic parameters derived from FAO, STATSGO, and SSURGO soil lookup tables, 
a neural network-based PTF (ROSETTA; Sehaap et al. 1998), PEST-Noah using PTFs, and 
measurements made during 2002 (Sehaap and Shouse, 2004) and 2004 (NAME). The FAO soil 
type for all 8 Metflux sites is sandy loam, STATSGO is loamy sand, and the finer resolution 
SSURGO data indicates 3 different soil types across the Metflux sites. As a result, there is 
significant disagreement in hydraulic properties amongst these lookup tables alone. 

The PEST-Noah parameters fall within the range of established datasets and measurements, 
yet there remain significant differences between lookup table and calibrated hydraulic properties. 
This is due to the ability of the PTF approach in Noah to result in a unique (but realistic) soil 
type that rigid lookup tables cannot describe. The ROSETTA PTF model suggests parameters 
that are inconsistent for a sandy soil and with observations, and indicate this particular PTF may 
not be appropriate for this region. Based on the improvements in simulated soil moisture 
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exhibited by PEST-Noah presented earlier, the PEST-Noah soil textures and PTFs appear to be 
the most accurate. 

It is also important to note that the large spread of hydraulic properties across data sources in 
Table 2 is a function of differences in way each property is estimated and what each represents. 
For example, K sat values of 250 cm d' 1 (NAME 2004) are estimated from soil samples under 
laboratory conditions, even though actual precipitation rates could never be high enough to 
observe similar saturation values in the field. Further, while the ROSETTA model suggests 
values for K sat that are an order of magnitude lower than those calculated in the laboratory (and 
may be closer to a true saturated value observed in nature), such values result in inaccurate 
simulations of soil moisture when employed in Noah. Overall, it is the combination of an 
accurate soil type representation that creates the range and scaling of hydraulic properties and the 
physics of the LSM that determines the most appropriate parameters in this ease. 

PEST-Noah suggests a unique soil type that also corresponds well with observations (Table 
2; Ksat and 8 S ). Although the model is not perfect, the limitations of the Noah soil physics are not 
significant enough to deter estimates of physically meaningful soil properties. On the contrary, 
previous attempts at model calibration in the region have yielded parameter values that lie well 
outside measured values, and as such could only be interpreted as ‘effective values’ that absorb 
significantly more inherent model or forcing data deficiencies than they do represent real soil 
properties. For example, Scott et al. (2000) found values of 2.5 x 10' 6 and 3.7 x 10 7 ms' 1 for K sat 
and 0.25 and 0.23 m 3 m 3 for porosity at Kendall and LH, using a model calibration approach. 
Using the PTFs, these K saf values correspond to a soil texture of 39 percent sand, and the porosity 
values are so low that a soil type cannot even be derived. 
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From a broader perspective, it is important to assess whether PEST -derived soil parameters 
can be employed in Noah and represent conditions at WGEW over longer timescales. Therefore, 
the soil textures optimized from the M9G period were used to run Noah at Kendall and LH over 
the 2002-4 evaluation period. Figure 6 shows the soil moisture simulated by PEST-Noah over a 
54-day period in the su mm er of 2003 compared with in-situ observations from Vitel probes 
surrounding LH and Kendall. Simulations with optimized parameters (RMSE = 2.4/2.1, Bias = - 
0.7/-1.9 percent for LH and Kendall) perform remarkably well compared to those using FAO 
soils (RMSE = 5.7/13.9, Bias = 5.6/13.3 percent) over the extended period encompassing 
numerous precipitation and drydown events. Once again, this highlights the ability of PEST to 
adjust the dynamic range of soil moisture simulated by Noah and effectively respond to 
precipitation events, and supports the use of optimized soil properties across this watershed for 
seasonal (and longer) durations. 

b. Temporal Sampling of PBMR Images 

The high temporal resolution of the PBMR imagery captures a complete soil drying cycle 
for this region. To assess the broader applicability of the methodology described above, it is 
useful to look at precisely how many and which PBMR images are needed for accurate 
calibration. The sensitivity of PEST-Noah to the number of PBMR observations was tested by 
testing at all possible combinations of image in the calibration process. Figure 6a shows the 
error in simulated versus observed soil moisture at Kendall for each of the image combinations 
used in PEST-Noah. Kendall was chosen as a representative site because it exhibits the largest 
range of soil moisture throughout the period and also was more difficult to calibrate due to 2 
significant precipitation events and a strong drydown in between. 
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These results show that there is a significant reduction in RMSE (and standard deviation) 
onee three or more images are used in the calibration. There are fewer number of data points for 
single images because many of those simulations were unable to converge with only one 
observation. What is also evident is there is a large amount of scatter or variability when using 
one or two images, but for three or greater all the points collapse indicating that it does not 
matter which images are included. Note that the error using Noah with FAO soils is over 2 times 
larger than even the worst PEST -Noah simulations using a single PBMR image. 

The other main factor in the success of PEST-Noah is what portion of the soil drying curve 
(i.e. dynamic range) is captured by the PBMR images. Figure 6b shows the error in simulated 
soil moisture against the range in soil moisture captured by the image combinations described 
above. The results look similar to Fig. 6a, and suggest that errors are significantly reduced if the 
images used capture at least 5 percent (volumetric) variability in soil moisture during a drydown 
period. When the foil dynamic range in soil moisture is captured by the PBMR images, the 
RMSE and bias (not shown) in PEST-Noah simulations are minimized, and are ~5 times lower 
than using FAO soils. 

Analyses also indicate that the second PBMR image (day 214) is the most critical 
observation to include in the calibration. This image was acquired immediately following a 
rainfall event (Fig. la) and represents the maximum value of soil moisture during the period. 
Out of the 7 simulations when PEST-Noah was unable to converge on a solution (i.e. not enough 
information was coming from the observations), all occurred when day 214 was not included. 
More importantly, out of the simulations using 5 out of 6 images, the only one unable to 
converge was with day 214 omitted. Error analyses (not shown) also support the importance of 
including day 214, and the improvement in calibrations when this ‘wet’ image is included. 
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An independent test of the sensitivity of PEST-Noah to the choice of PBMR images was also 
conducted. PEST-Noah calibrations were performed for the early part of the M9G period using 
only the first 2 or 3 PBMR images, then evaluating Noah over the remainder of the period with 
the optimized soil texture as input. The RMSE in simulated soil moisture using this approach are 
plotted alongside the dependent results in Fig. 7, and are below 2 percent and similar to those 
from calibrating and evaluating the whole M90 period. This suggests that on very short temporal 
scales, PEST-Noah could be used with a few images early in the period to calibrate and estimate 
soil texture, which could then be used to improve simulations going forward without requiring 
additional images. This type of approach mimics that of an operational and data-limited 
approach such as that of the ARMS project. 

Overall, these results demonstrate that the majority of the improvement using PEST comes 
from calibrating the overall magnitude and dynamic range of soil moisture by adjusting the sand 
and clay contents. Also, PEST-Noah significantly reduces errors in simulated soil moisture 
regardless of the number and which images used. However, the calibration is most accurate 
when you include more than two images in the calibration process, or at least one or two images 
that capture and represent the observed range of soil moisture during a drying cycle. 

c. 2003 Calibration Experiments 

The development of the delta index allows us to test the PEST-Noah approach using active 
radar imagery. RADARSAT-1 images were acquired over the WGEW during July, August, and 
September 2003 that cover a larger temporal and spatial extent than the PBMR images during the 
M90 period (Fig. 8). As described earlier, the Delta Index data had to be aggregated from 40 to 
210 and 280 meter resolution to reduce the impact of speckle (amplified by the high rock content 
of the soils in this region) on the soil moisture retrievals. PEST-Noah was still run with 
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precipitation and land cover data at 40 meter resolution, and therefore required experiments to be 
conducted to confirm the appropriateness of applying 210 and 280 meter resolution data from 
active radar to that of a 40 meter pixel. 

Analyses showed that the maximum differences in observed PBMR soil moisture within a 
280 m 2 pixel surrounding each Metflux site were less than 2 percent and well within instrument 
error (± 5 percent). PEST-Noah was then run individually over 7x7 points (to match the 
resolution of the delta index) surrounding the Kendall site during the M90 period using the 40 
meter PBMR data to calibrate. Even at the site with the largest precipitation, land cover, and soil 
moisture gradients (Kendall), the degree of spatial heterogeneity on the scale of 280 m 2 does not 
lead to significant differences in soil properties estimated by PEST-Noah. The remaining sites 
are even more uniform in nature, which therefore gives confidence to using the 210 and 280 
meter delta index forLH and other sites in 40 meter PEST-Noah calibrations. 

To cover the extended period between active radar overpasses, PEST-Noah was run from 30 
June - 15 September 2003 using the interpolated precipitation and Kendall forcing data as for the 
2003 evaluations described previously. For the Kendall site, PEST-Noah converged on a 
solution of 100 percent sandy soil for both the 210 and 280 meter images (compared with 89 
percent sand; 1 1 percent clay in M90). At LH, the optimized values for sand, clay, and silt were 
28, 72, and 0 percent using 210 meter data, and 20, 45, and 35 percent using 280 meter data 
(compared with 99 percent sand in M90). 

Figure 9a shows that results from Kendall are comparable to the in-situ soil moisture 
observations, and although this calibration suggests a slightly sandier soil than M90, the 
differences in soil moisture and hydraulic properties are negligible as discussed in Section 3a. 
Following the discussion in Section 3b, the differences in the M90 and 2003 period calibrations 
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are primarily due to the lower magnitude range of soil moisture captured by the delta index 
(0.153 max, 0.032 min) compared to the PBMR (0.169 max, 0.075 min). That the first image 
coincides with a rainfall event means the soil type is less relevant or sensitive to calibration on 
this date, and actually reflects the need to dry out to match the very low soil moistures given in 
the next two images. 

The PEST-Noah results for LH (Fig. 9b) suggest considerable differences in ealibrations- 
using the delta index compared with observations and the PBMR calibrations, with a much more 
clayey soil and higher soil moisture throughout the period. Closer inspection of the in-situ 
observations of soil moisture over the period shows a comparable range to Kendall. However, 
the first delta index image (at both 210 and 280 meters) gives a rather wet soil moisture estimate 
compared to observations, while the latter 2 images are relatively dry and alone would suggest a„ 
sandy soil as in previous calibrations. In order to match the high moisture content of the first 
image, PEST-Noah is forced to simulate a high clay content which has a higher holding capacity 
and strength. 

The first radar overpass is more than 4 days after the last significant rainfall, which means 
the soil has had significant time to dry out particularly for this region. In fact, studies have 
shown that the typical response time to rainfall and complete drydown occurs within 2-3 days (or 
less, depending on ground cover) in WGEW due to the shallow moisture reservoir and high bare 
soil evaporation rates (Shamir et al. 2005; Kure and Small 2004). That previous studies, M90 
observations (PBMR and in-situ), and in-situ observations all depict a much more rapidly drying 
soil at the site suggests that the delta index data may not be accurate on this date. 

High frequency (active) microwave retrievals are difficult to obtain in regions with high 
rock content due to increased backscatter and a weaker relationship with soil moisture (Jackson 
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1992) even when aggregated to 280 meter resolution, and thus remains an issue. Clearly, if there 
were more images available, particularly during and immediately following rainfall events, 
PEST-Noah would be able to perform better as for the M90 case. Houser et al. (1998) made 
similar recommendations for data assimilation in this region, suggesting that soil moisture 
observations are required at minimum once every storm event. For the 2003 experiments, it is 
likely a combination of insufficient temporal sampling and the limited spatial resolution of active 
remote sensing (through the signal-to-noise ratio) in the soil moisture retrieval process that 
resulted in poor calibration at LH. 

5. Discussion 

As detailed point or regional surface characteristics are not always available, it is useful to 
examine the calibration technique at lower spatial resolution. Figure 10 shows simulated soil 
moisture from default (FAO) Noah and PEST-Noah simulations compared with the PBMR data 
for the entire PBMR domain on DOY 221. In effect, PEST-Noah is minimizing the mean error 
in simulated versus observed soil moisture across this region. The last PBMR date was chosen 
because it represents the cumulative effect of calibrated parameters on improving soil moisture 
throughout the M90 period, lowering both RMSE and bias from ~10 to 3 percent. While there is 
only one set of optimized soil texture and hydraulic parameters estimated for the entire domain 
(92 percent sand), the significant improvements indicate that on watershed scales with coarser 
inputs the calibration process can still successful and potentially useful for deriving meaningful 
soils data. 

Overall, PEST-Noah is able to identify the dominating soil texture and hydraulic properties 
of the WGEW. It is also important to emphasize that the calibration approach was designed to 
focus on calibrating a consistent set of hydraulic parameters that would retain the physical 
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characteristics of the region in question. Sophisticated techniques have been successful in 
calibrating complex models, but yield ‘effective’ parameters that have lost physical meaning. 
This practice is troubling for the future development and utilization of complex LSMs with 
respect to the applicability of new and improved parameter data that are likely to be available 
from remote sensing techniques. As LSMs become increasingly complex, the ability to retain 
physically meaningful and measurable parameters needs to be addressed in parallel, in a manner 
similar to this study. 

However, it is accepted that some of the errors in model physics, forcing, and 
parameterizations are accounted for in PEST-Noah calibrations as well, and should be addressed 
(Scott et al. 2000). While the approach taken here is rather simple compared to complex data 
assimilation or multi -objective parameter estimation techniques, it is easier to identify errors in 
the model. One deficiency in the Noah model was identified, where the bare soil evaporation 
factor should be adjusted (lowered) to allow for greater evaporation rates and soil drying for 
semi-arid regions with sparse vegetation. The original value, when combined with FAO, 
STATSGO, and SSURGO parameters in the default simulations led to soils that were 
consistently too wet for this region (as shown in the results). 

It has been noted in numerous studies of the WGEW that there is an unusually high rock 
fragment content of the upper soil layers (and is not typically accounted for in categorical lookup 
tables or PTFs). The unique soil type estimated by PEST-Noah at the Metflux sites is a primarily 
sandy soil that is slightly less porous and more conductive than that observed. In fact, this may 
be an attempt to indirectly account for the high rock content in Noah by adjusting the parameters 
to match that of a rocky, yet sandy, soil (i.e. less pore space and increasing flow paths in the soil 
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volume). A simple formulation to account for rock content in calculating hydraulic parameters is 
currently being implemented into Noah to test the sensitivity of the calibrations. 

These issues demonstrates the potential for remote sensing data to also offer information on 
model errors and biases (although small in this case), and advanced techniques such as those 
offered by Bach and Mauser (2003) and Ines and Droogers (2002) to find and quantify 
systematic errors in LSM physics. In particular, the recent and ongoing work of Amramowitz et 
al. (2006) and the Model Parameter Estimation Experiment (MOPEX; Duan et al. 2006) are 
promising new avenues of research that are working to determine the applicability of studies 
such as this one to be applied to other LSMs and locations. 

6. Conclusions 

This paper has examined a straightforward method of using microwave remote sensing of 
near-surface soil moisture to calibrate an offline land surface model, and in the process infer soil 
texture and hydraulic properties at high spatial resolutions. This approach expands and improves 
upon a wide body of previous work by incorporating pedotransfer functions into the LSM to 
ensure consistent and physically meaningful soil parameters, and by addressing the temporal 
sampling of remote sensing imagery needed for successful calibration. As a testbed for the 
AJRMS project, this research was able to retrieve soil texture and property estimates that 
correspond well with observed soils over the WGEW. Once estimated for this region, these 
parameters were also used to simulate soil moisture over seasonal time scales with a great deal of 
accuracy compared to simulations with default soils and soil properties based on lookup tables. 

Specific results of this study include the following: 

1) Limited microwave retrievals of near-surface soil moisture can be used to calibrate a 
LSM to within .02 m 3 rrf 3 accuracy at high temporal and spatial resolutions. 
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2) Optimizing soil hydraulic properties using PTFs gives better and more physically 
meaningful results than a one-at-a-ti me parameter estimation approach. 

3) Errors in the calibration process are minimized when there are at least 3 images included 
that represent the typical range of moisture exhibited by the soil type during a drydown 
period. 

4) Independent tests indicate that this methodlogy can be successful in calibrating LSMs 
over seasonal and longer timescales for use in specialized prediction systems. 

Overall, these results suggest that ARMS could be applied at remote locations to simulate 
soil moisture in a semi-operational context with limited remote sensing inputs. Simulations that 
expand the 8 Metflux sites tested here to the full WGEW at 40m resolution of PEST-Noah are 
ongoing, from which fully distributed maps of soil texture and hydraulic properties will be 
produced. Alternatively, one could use the PEST-Noah approach after stratifying the watershed 
using high resolution soils, land cover, or similar data to further examine the spatial distribution 
of soil properties. Distributed soil property information can then be compared with existing soils 
maps and the approach repeated and applied to other LSMs and regions. 

Finally, the ability of active remote sensing and the delta index technique to retrieve soil 
moisture on less than watershed scales needs to be investigated further before incorporated in an 
ARMS-type of approach. While results here have shown that 3 images are sufficient to calibrate 
and obtain soils information, the soil moisture estimates must be accurate (within ~3-5 percent) 
and capture a typical dynamic range of soil moisture for the region in question. The spatial 
resolution of currently orbiting active remote sensing, determined in part by the signal-to-noise 
ratio of the measurement, may be a limiting issue for this application. Also, the accuracy of soil 
moisture retrieval from active remote sensing through approaches like the delta index or other 
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retrieval methods (e.g. Alvarez-Mozos et al., 2005) needs further investigation before it can be 
incorporated with confidence. Future work on the ARMS project will include testing the 
methodology and evaluating the delta index at cold land, high relief, and highly coupled regions 
of the U.S., where the calibration process should yield new insight [about?] the images and 
accuracy required and influence of model physics for diverse surface conditions. 
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Table 1: Optimized sand, silt, and clay percentages estimated from a) PEST-Noah 
simulations at the eight Metflux sites compared with b) those observed by Schmugge et 
al. (1994) and during the NAME (Higgins et al. 2006) in 2004, and associated hydraulic 
properties computed for each using the PTFs employed in Noah (Cosby et al. 1984). 
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percentages at each site using the PTFs in Noah. For comparison, soil properties estimated by Scott et al. (2000) from a neural 
network PTF (ROSETTA; Schaap et al. 1998) are shown along with site-specific estimates of hydraulic parameters based on soil 
samples taken during 2002 (Schaap and Shouse 2004) and NAME (Higgins et al. 2006) in 2004. 
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Figure 1: The Walnut Gulch Experimental Watershed in southeastern Arizona (outlined 
in black) covers 148 km 2 and is heavily instrumented with meteorological, flux, and rain 
gauge data. The M90 experiment included 8 Metflux sites (°) of which Lucky Hills (Site 
1) and Kendall (Site 5) were supersites. Overlain are estimates of volumetric soil 
moisture (m 3 m' 3 * 100) derived from Push Broom Microwave Radiometer measurements 
onDOY 214. 
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Figure 2: Soil moisture in the 0-5 cm layer at the a) Kendall and b) Lucky Hills sites 
during the Monsoon ’90 study period from (■) PBMR retrievals (Schmugge et al. 1994) 
and (A) gravimetric measurements with standard deviations of the 3 measurements made 
at each site. Also plotted are the 6-hourly precipitation totals during the period at each 
site as derived from the 84-gauge interpolated dataset over WGEW (Garcia et al. 2006) 




Site 5 - Kendall 


B 

PBMR 

A 

GRAV 

- - 

■FAO 

L r . J 


i' 


lf)<DN-C0050T-OJ« r 5’<rU2^N : OO^OT^ 

ojoJOJ^Soj<NC'l<NOJC'JC'i<NC'J040tOJ 


Site 1 - Lucky Hlls 


-a — PBMR 
A GRAV 

FAO - 

-» — PEST 


iocoN-coa>o*-o4cO'sj-ioc0N-coa>o 

O0000^-^“T"^” T-T-T “ T " T_T ~ C N 

<N<N<N<NCs1C\JC\1<N<NO'1CN<NCNC9CNCM 

Figure 3: Simulated 0-5 cm soil moisture from default and PEST-calibrated Noah 
simulations for the a) Kendall and b) Lucky Hills sites during the M90 period. 
Measurements of soil moisture from PBMR and gravimetric sensors on the 6 
optimization dates are also shown. 
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Figure 4: Bias and RMSE in simulated versus observed (PBMR) 0-5 cm soil moisture 
during the M90 period using default (FAO; gray) and optimized (PEST; black) soil 
properties at each Metflux site. 




Figure 5: Percentages of sand, silt, and clay estimated by PEST-Noah at the eight 
Metflux sites compared with in-situ measurements from Schmugge et al. (1994). 
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Figure 7: Errors in simulated versus observed 0-5 cm soil moisture at the Kendall site for 
varying a) numbers of PBMR images used in PEST-Noah and b) ranges of soil moisture 
covered by these images. The lightly shaded points indicate simulations that were 
calibrated using only the first two and three PBMR images. 
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Figure 8: Soil moisture (m3 m 3 ) estimated from RADARSAT-1 active microwave 
measurements over the WGEW on a) 29 July, b) 22 August, and c) 15 September 2003. 
Backscatter was aggregated from 7 to 280 meters to reduce the effects of speckle and 
converted to soil moisture using the delta index image differencing technique (Thoma et 
al. 2006). 
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Figure 9: Near-surface soil moisture simulated by PEST -Noah using 210 and 280 meter 
Delta Index as observations at the a) Kendall and b) Lucky Hills sites during summer 
2003. Also shown are observations of soil moisture from Vitel probes surrounding each 
location. 






